The parametric oscillation threshold of semiconductor microcavities in the strong 
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The threshold of triply resonant optical parametric oscillation in a semiconductor microcavity 
in the strong coupling regime is investigated. Because of the third-order nature of the excitonic 
nonlinearity, a variety of different behaviours is observed thanks to the interplay of parametric 
oscillation and optical bistability effects. The behaviour of the signal amplitude and of the quantum 
fluctuations in approaching the threshold has been characterized as a function of the pump, signal 
and idler frequencies. 
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I. INTRODUCTION 

Triply resonant optical parametric oscillation 
(OPO)ii£ has been recently observed 3 -!^ in semi- 
conductor microcavities in the strong coupling regime^ 7 -, 
and has attracted a good deal of attention from the 
point of view of both fundamental physics and possible 
technological applications. The peculiar dispersion 
relation of polaritons in the strong coupling regime 
allows to simultaneously satisfy the resonance condition 
for the pump, signal and the idler modes. Together 
with the enormous value of the excitonic nonlinearities, 
the possibility of easy phase matching results in a low 
threshold intensity, making these systems very promising 
candidates for low-power OPO applications. 

A complete theoretical description of the OPO dynam- 
ics of such systems is not only very important in view of 
the optimisation of the device operation, but also de- 
serves a certain interest from the point of view of nonlin- 
ear dynamics as many interesting phenomena can occur 
due to the interplay of optical bistability and parametric 
oscillation 8 -^, and to the nontrivial spatial field dynamics 
in the transverse pian o 10 ' 11 . 

As shown by several theoretical papers that have ap- 
peared on the subject, a rather complex phenomenology 
is found already at the level of the three-mode approxi- 
mation, where the classical nonlinear optical wave equa- 
tion is projected onto the three pump, signal and idler 
modes^ ' 12 ' 13 ' 14 . Available experimental data appear to 
confirm this point: in particular, both continuous^^ 
and discontinuous^ 3 - behaviours have been experimen- 
tally shown for the signal intensity in the neighborhood of 
the threshold point. Although some analogies have been 
drawn with what is known about x OPO dynamics 
in standard passive medial ' 16 ' 17 ' 18 , no complete inves- 
tigation has appeared yet for the case of semiconductor 
microcavities in the strong coupling regime, neither from 
the experimental nor from the theoretical points of view. 

The optical nonlinearity of the microcavity system 
under investigation originates from collisional exciton- 
exciton interactions and is therefore of the x tyP e - 
This means that it not only provides the parametric in- 



teraction necessary for the parametric oscillation, but is 
also responsible for significant mean-field frequency shifts 
of the modes. This makes the nonlinear dynamics of the 
mode amplitudes much richer than in x^ 1 OPOsi 8 -. Pi- 
oneering theoretical work in this direction has recently 
appeared in Ref. 0. 

The purpose of the present paper is to provide a sys- 
tematic and quantitative study of the OPO threshold 
in semiconductor microcavities in the strong coupling 
regime. Depending on the pump laser frequency, and the 
signal, pump and idler mode frequencies, several regimes 
are to be distinguished, where the system behaviour is 
radically different. 

The paper is organized as follows: our model of the 
microcavity is introduced in SecJTTl Optical limiting and 
optical bistability in the pump only solution are discussed 
in Sec lIII Al General concepts about the stability of the 
solution with respect to pump-only and to parametric in- 
stabilities are given in Sec lIII Bl The following Sccs lIII CT 
I1II Gl are devoted to characterize the parametric thresh- 
old as a function of the incident pump angle, the internal 
and the incident intensities and to find the optimal choice 
to minimize the threshold intensity. Quantitative esti- 
mations are provided in Sec lIII Hi where a comparison 
is made with other realizations of OPOs based on pas- 
sive x^ an d x'" 3 '' materials. The kind of bifurcation at 
the onset of the OPO emission is investigated in Sec lIVI 
Depending on whether the pump-only solution is in the 
optical limiter or in the optical bistability regimes, para- 
metric emission is shown to set in either in a continuous 
or in a discontinuous way. The close relationship between 
the nature of the instability point and the behaviour of 
the quantum fluctuations as the threshold is approached 
is pointed out in Sec lIVDl Conclusions are finally drawn 
in SecEl 



II. POLARITON MODEL 

A sketch of the physical system under investigation 
is shown in FigfT] a planar DBR (Distributed Bragg 
Reflector) semiconductor microcavity containing a few 
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quantum wells strongly coupled to the cavity mode. The 
elementary excitations of this system consist of exciton- 
polaritons, i.e. coherent superpositions of cavity photons 
and excitons. Polaritons combine the very strong x*- 3 -* 
optical nonlincarity originating from exciton-exciton col- 
lisional interactions to the peculiar dispersion relation as 
a function of the in-plane wavevector k that is shown 
in FigfT] these facts make them extremely well suited 
for triply-resonant optical parametric oscillator applica- 
tions, as it has indeed been experimentally demonstrated 
in recent years 3 -*^. 
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FIG. 1: Upper panel: sketch of the microcavity system and 
of the parametric process under consideration. Lower panel: 
lower polariton (LP) and upper polariton (UP) dispersion 
at linear regime. The cavity photon dispersion is ujc{k) = 
ujc ^l + k 2 /k 2 z with Tiu)% = 1.4 eV and k z = 20 /an -1 . The 
exciton dispersion is flat and resonant with the k = cav- 
ity mode lux = Uc- The exciton- photon Rabi coupling is 
hQ.R — 2.5 meV. The dots indicate the signal, pump and idler 
modes, the arrows show the triply resonant parametric pro- 
cess under investigation. 



A mean-field description of the cavity-polariton field 
dynamics can be developed in terms of a nonlinear wave 
equation with a third-order nonlinearity^i 3 -. Under the 
assumption that the dynamics takes place in the lower 
polariton branch and the population of the upper polari- 
ton branch remains negligible, the theoretical description 
can be simplified by restricting it to the lower polariton 
only. For the sake of simplicity, we shall focus our at- 
tention on the case of a circularly polarized pump beam. 
As the circular polarization of the polariton field is pre- 
served by the nonlinear interactions and the longitudinal- 
transverse splitting 19,20 is much smaller than both the 
linewidth 7 and the nonlinear interaction energy, the cir- 
cular polarization is almost completely transferred to the 



signal and idler beams. 

Under these assumptions, the polariton dynamics can 
be written in terms of a nonlinear wave equation for a 
single-component fc-space polariton field ^/>£p(k): 



e{k) - i 



7 (k) 



^LP(k)+F p (k) e 



qi,q2 



(1) 



The field V'Lp(k) is here normalized in such a way that 
its square modulus |V'Lp(k)| 2 equals the number of po- 
laritons with momentum k per unit area. e(k) is the 
dispersion relation of the lower polariton and v(k) is the 
momentum-dependent loss rate. Throughout the present 
paper, the exciting laser field is taken as a monochro- 
matic and continuous wave coherent field at uj p with a 
plane-wave spatial profile at k p and a circular polariza- 
tion. The driving amplitude F p (\s) can be related to the 
incident power density 7j„ c by using the input-output 
formalis m 21 ' 22 i 23 : 



F p (k) = <5k,k p C(k p 



/ ^radXim 
N+rhhJv 



(2) 



^rad is here the radiative decay rate of the cavity-photon 
due to the finite mirror transmittivity; the parameter 
Ntr specifies whether the cavity is a single-sided cavity 
with a perfectly reflecting back mirror (Ntr = l)j or a 
symmetric cavity with equal transmission through both 
the front and back mirrors (N tr — 2). 

The third-order nonlinear interaction term takes into 
account exciton-exciton collisional interactions. As the 
wavevectors involved in the present discussion are much 
smaller than the inverse excitonic radius, the exciton- 
exciton coupling constant in a single quantum well can be 
approximated by a momentum- independent g. If Nqw 
quantum wells are present in the cavity, identically cou- 
pled to the cavity mode, the bright excitonic excitation 
is delocalized over all of them and the effective excitonic 
coupling constant is g = cj/Nqw- In the polaritonic ba- 
sis, a non-trivial momentum dependence appears via the 
Hopfield coefficients A(k) and C(k) quantifying the ex- 
citonic and cavity-photonic components of the lower po- 
lariton: 

ff k , qi , q2 =ffX*(k)X*(q 1 + q 2 -k)X(q 1 )X(q 2 ). (3) 

Although no conclusive experimental nor theoretical 
analysis has been reported yet, the theoretical predic- 
tion Hg s» 1.5 x 1CU 5 eV fim 2 based on the Born approx- 
imation 24 appears to be in reasonable agreement with 
available experimental data^i. 

In order to focus our discussion of the basic OPO dy- 
namics, we shall not consider here the effect of the disor- 
der: in recent high quality III-V samples the effect of the 
disorder can in fact be weak enough for it to be neglected 
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on the scale of the polaritonic linewidth. In this case, it 
is legitimate to approximate the mode eigenfunctions as 
plane waves. On the other hand, the disorder is much 
stronger in II- VI samples, where it has been shown to 
have dramatic consequences on polariton BEC— . These 
effects are highly non-trivial already at equilibrium^ 6 , and 
are expected to become even more complex because of 
the interplay with the nonlinear dynamics: the complete 
analysis of them goes far beyond the scope of the present 
paper and is left to future work. 

To conclude the section, it is interesting to note that 
the applicability range of the wave equation (fT]) is not 
limited to semiconductor planar OPOs, but can be ex- 
tended to describe other setups, e.g. planar cavities con- 
taining a slab of passive x^ material. In this case, no 
excitonic resonance exists, and the polariton reduces to a 
bare cavity-photon. As both the coupling to external ra- 
diation and the optical nonlinearity act on the same pho- 
tonic degree of freedom, one has simply to set X = C = 1 
and calculate the nonlinear coupling constant using the 
nonlinear susceptibility of the medium under considera- 
tion: 

K 9 = C X W^. (4) 

The numerical factor C of order one takes into account 
the details of the geometry under investigation. Typical 
values of x of materials specifically designed for non- 
linear optical applications range up to something of the 
order of 10 _9 esv* 2 ^.. For a A/2 cavity, these values cor- 
respond to a nonlinear coupling constant of the order of 
hg w 5 x 10 _9 eV/xm 2 , orders of magnitude lower than 
the value hg ss 1.5 x 10 -5 eV/im 2 previously mentioned 
for semiconductor microcavities in the strong coupling 
regime. This explains the present interest of semiconduc- 
tor microcavities for low-power nonlinear optical applica- 
tions, as well as for fundamental studies of the interplay 
of nonlinear dynamics and quantum fluctuations^^ 9 -. 

III. THE PARAMETRIC THRESHOLD 
A. The pump only solution 

Among the different k modes, only the one at k = k p 
contains a source term in its equation of motion (TT]). An 
exact solution of the full set of equations of motion (fT|) 
can therefore be found in the form 

fe^il^Pe'^'', (5) 

with the amplitude P fixed by the condition: 

[e p - wp - % - lp + 5 |X(fc p )| 4 |P| 2 ]P + F p = 0. (6) 

Here e p = e(k p ) is the frequency of the pump mode 
at linear regime and j p = 7(k p ) the corresponding 
linewidth; the effect of the third-order nonlinearity is 
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FIG. 2: Polariton density in the pump mode as a function 
of the driving intensity I p = |P P | 2 for a pump wave vector of 
k p = 1.19/im -1 , which corresponds to e p = 1.39845 eV. The 
left graph is in the optical limiter regime (u> p = e p ), whereas 
the right one is in the bistability regime (ui p = e p + 1.57 P ). 
The red arrows show the jump in polariton density for an 
upward (right arrow) and downward (left arrow) ramp of the 
laser intensity. The dotted part of the curve is dynamically 
unstable. Same cavity parameters as in Fig[TJ 



to renormalize the pump mode frequency by a mean- 
field shift proportional to the mode excitonic population 
n xp = |A(k p )| 2 |P| 2 . This effect, absent in x^ cavi- 
ties, is responsible for the the qualitatively different be- 
haviours^ that can be observed depending on the sign 
of the detuning of the pump frequency lu p with respect 
to the polariton energy e(k p ). 

FigO shows the excitonic density in the pump mode 
function of the driving intensity I p = \F p \ 2 which 
is proportional (but not identical) to the incident laser 
intensity i,- nc [See Eq. (|2j)] . When the pump frequency is 
below or close to resonance u> p — e(k p ) < \/37(kp)/2, we 
are in the so-called optical limiter regime, in which the 
population n xp of the k p mode monotonically increases 
as a function of the driving intensity I p (left panel). 

For blue-detuned pump frequencies lo p — e(k p ) > 
V / 3 7(k p )/2, a positive feedback of the nonlinearity oc- 
curs and hysteretic behaviour can be instead observed, 
as shown in the right panel of Fig[2] and experimentally 
demonstrated in Ref. l3ll . For increasing laser intensity, 
the pump mode population follows the lower branch until 
its endpoint is reached, and then it jumps to the upper 
branch as indicated by the arrow. If the driving inten- 
sity is later ramped down, the system keeps following the 
upper branch until its endpoint, and only here the pump 
mode population jumps back to the lower branch. Hys- 
teretic behaviour is apparent, as the upward and down- 
ward jump points do not coincide. 

B. Dynamical stability of the pump-only state 

As usual in nonlinear dynamical systems, finding a so- 
lution is not sufficient, as one has to verify its dynamical 
stability. Optical parametric oscillation, as well as the 
instability of the central branch of the hysteresis loop 
are in fact due to the solution ^ becoming dynamically 
unstable. 
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As the planar cavity supports a continuum of inde- 
pendent modes with different in-plane wavevectors, dy- 
namical stability of the pump only solution ([6]) has to 
be checked with respect to perturbations with any wave 
vector k s : 

7/>Lp(x, t) = p e -^pt+*px + u ^ e -i[w p +w(k a )]t+ik e x 

+ t!*(k s )e _4 '" p ~^*^ ks ^ t+ ^ 2kp_k ^ x . (7) 



Substituting this expression in Eq.([T]) and keeping only 
linear terms in the fluctuations u and v, one gets to the 
following eigenvalue problem^ 

L(k s )w(k s ) = u(k s )w(k s ) (8) 

where the two-component vector w(k s ) = [u(k s ), v(k s )] T 
and the 2x2 matrix L(k s ) is 



T(v ^_ (e s -u p + 2 9 \x s \ 2 \x p \ 2 \p\ 2 -if gx;XiX*i* \ 

L{ks) { -gX s X*XfP^ -e t + us p - 2g\X l \ 2 \X p m 2 -if)' (9) 

The matrix L(k s ) couples the fluctuations in the k s and k^ = 2k p — k s modes, called in the following the signal and 
the idler modes. Short-hand notations have been here introduced to simplify the expressions: X s ,i = X(k Sj i) are the 
excitonic Hopfield coefficients of the signal/idler modes, e S: i — e(k S) j) are the signal and idler mode frequencies and 
7s, i = 7(k s ,i) are the corresponding loss rates. Dynamical stability is ensured if the imaginary parts of all eigenvalues 
of L(k s ) are negative Im[w±(k s )] < for all wavevectors k s . These can be written as: 




Im{c±(k s )} = ^ll±Jl ±Im J [(e si - w p ) + g(\X s \* + \X^)n xp - i^-^L]2 - g * \Xtf\Xtfn* , (10) 



Note that the pump mode frequency e p does not directly 
appear in the expression (|10[) of the eigenvalues, but it is 
only indirectly involved via the pump-only solution 
which fixes n xp . The frequencies of the signal and idler 
modes are involved in (|10[) only via their average value 
(-si = [e a + eJ/2. 

Two kinds of physically distinct instabilities can arise. 
A single-mode instability arises when the equation of mo- 
tion for the pump mode alone - neglecting all interac- 
tions with other modes - is dynamically unstable. This 
instability is found when L(k s ) has an eigenvalue with 
a positive imaginary part for k s = k p . As in this case 
hi = k s = k p , this instability involves the k p mode only, 
and for this reason it is called single-mode. It is easy 
to verify^ that the pump-only solution © is single-mode 
unstable in the central branch of the bistability curve 
(marked with a dotted line in FigJ^b). At the turning 
points of the bistability curve a stable and an unstable 
solution meet, so that the bifurcation is of the saddle 
node type^. 

Our interest is however more focussed on instabilities 
of the second kind, i.e. for k s ^ k p : this parametric in- 
stability signals the onset of parametric oscillation with a 
finite intensity appearing in a pair of distinct signal/idler 
modes at k Sj i. From the point of view of bifurcation 
theory, the parametric instability profoundly differs from 
the single-mode one. As we shall see in the following, 
the pump-only solution still exists beyond the threshold 
point, but it is no longer stable for an eigenvalue of the 
linear stability matrix L(k) has crossed the real axis: the 
bifurcation is then of the Hopf type 33 and is accompanied 



by a spontaneous breaking of a signal/idler U(l) phase 
rotation symmetry^. 



C. Available range of signal/idler frequencies 

In the present paper, we shall not address the problem 
of the determination of the wavevectors k Sj i which are 
actually selected by the parametric process above thresh- 
old. This is a very complicate problem and is postponed 
to a forthcoming publication 11 . Here we shall limit our- 
self to a study of the lower threshold for parametric emis- 
sion: the parametric oscillation dynamics will be initiated 
as soon as the incident intensity exceeds the threshold 
value for some pair of signal/idler modes. 

For each value of pump wavevector k p , it is important 
to characterize the range of e S i that can be obtained when 
the signal/idler wavevectors k Sj i are spanned through all 
different polariton states: the search for the minimum 
value of the threshold has in fact to be restricted to the 
region of e s , values which are actually available. 

This point is addressed in Fig[3J In the left panel, 
the behaviour of the detuning (e s j — e p ) as a function 
of k s is shown for three different values of k p and the 
yellow region in the right panel summarizes the accessible 
detunings as a function of k p . For small k p , the e S i vs. k s 
curve has a single minimum at k s = k p where e S i = e p , 
and then tends to a finite limit for large k s (e(fc) has a 
finite limit for large k). For larger values of k p , negative 
values e s i — e p < can be reached. The minimum is in 
fact split in two separate minima^, symmetrically located 
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FIG. 3: Left panel: plot of e s ; as a function of k a for fixed 
values of k p = (full black), k p — 1.19/^m - (dashed red), 
k p — 1.4 fxm~ (dotted blue). The vertical lines show the 
value of the pump wave vector. Right panel: band of available 
e s i values as a function of k p . The horizontal green line shows 
the optimal detuning A°f^ (see Sec lIIIFl below), the vertical 
lines indicate the k p values corresponding to the curves in the 
left panel. Cavity parameters as in FigfTJ equal damping rates 
"fs. P ,i = 7 with tvy = 0.25 meV. 



around the pump angle as required by the symmetry of 
e S i under exchange of the signal and idler modes. The 
upper limit of the available band monotonically decreases 
as a function of k p , due to the corresponding increase of 
e p . In particular, it tends to for large values of k p . 
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FIG. 4: Threshold excitonic density in the pump mode as 
a function of the pump laser detuning from the average sig- 
nal/idler natural frequencies. The full black curve refers to 
a semiconductor microcavity for which the Hopfield coeffi- 
cients are X s = X(0) and Xi — X(2fc mag i c ). The dashed 
red curve refers to the case of a cavity containing a passive 
nonlinear material for which Xj = 1. The curves do not con- 
tinue through the hatched region at the left hand side where 
parametric oscillation can never take place. The dotted lines 
represent the approximation (|14p . Cavity parameters as in 
FigH 



D. Pump intensity n xp at the parametric threshold 

As it often happens in nonlinear optical systems, it is 
useful to study the parametric threshold first in terms 
of the internal light intensity in the cavity, in our case 
the excitonic pump mode population n xp . Connection to 
the incident intensity Ii nc will be then made in the next 
subsection. As we are still left with several parameters, 
namely Ji/f s , X 2 anci Xf, we are forced to restrict the 
discussion to some illustrative examples. The qualitative 
features are however quite robust with respect to their 
variations. Let us begin from the 7 S = 7i = 7si case: as 
the argument of the square root in eq. (|10p is purely real, 
the calculations are in this case the simplest. 

The pump mode population n xp at the parametric 
threshold is plotted in Figf4] as a function of lo p — e S i 
for two possible choices of the Hopfield coefficients. No 
qualitative differences are visible, but only quantitative 
ones. The main feature of these curves is the fact that 
parametric oscillation can only take place for sufficiently 
large values of cu p — e S! ;. The hatched regions indeed in- 
dicate where parametric oscillation can never take place, 
no matter how large the population of the pump mode is. 
Remarkably, the minimum value of the threshold popu- 
lation is reached just before the endpoint of the curves. 

Simple physical arguments can be put forward to ex- 
plain these features. In a v/ 3 -* parametric oscillator, the 
nonlinearity not only provides the parametric coupling 



between the signal and idler modes via the off-diagonal 
terms in the matrix ([9]), but is at the same time responsi- 
ble for a mean-field blue shift of the signal and idler mode 



frequencies by 2 g \X. 



s,i\ ,L xp 



Once this shift is taken into 
account, the resonance condition for the parametric pro- 
cess is renormalized to 



top = £ S i + g{\X s \ 2 + \Xi\ 2 )n xp . 



(11) 



From (|10p . it is easy to see that the minimum value of 
the threshold 



7s 



Xf) 



2g\X s \ \Xi 



(12) 



is indeed attained when this condition is satisfied. Com- 
bining (jlip and (TT2"]) gives the optimal detuning 



IX., 



\Xi 



7s 



(13) 



which corresponds to the position of the minimum of the 
curves plotted in FigdJ 

On the other hand, for large and positive values of the 
detuning lo p — e S i 3> 7 S i, the threshold grows in a linear 
way as a function of lo p — e S i 



l xp 



gdX^ + lX^ + lXsX^)- 



(14) 



Finally, for lj p — e S i < the well-known inequality 
(|A S | 2 -I- \X,\ 2 ) 2 > \X S \ 2 IX,] 2 implies that dTDJ) can never 
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FIG. 5: Contour plot of the threshold laser intensity as a 
function of the detunings. The decay rates for the pump, 
signal and idler modes are taken equal to 7 s , P ,i = 7 and 
the Hopfield coefficients X a = X(0), X p = X(fc mag ic) and 
Xi = X(2fc mag i c ). The lowest contour line in the plot is 
at 1.05/j™™, where the minimum of the threshold intensity 
1™™ is attained at the point indicated by a star * and is de- 
fined in Eq. (|f 8|l below. The difference between the contours 
is 0.3/j™"- The letters (A)-(E) indicate the regions of qual- 
itatively different behaviours; the corresponding pump-only 
characteristic curves are shown in Fig(S] The red-shaded area 
indicates the values of e s ; — e p that are available for the value 
of the pump angle k p — 1.4pm , corresponding to the red 
curve in Fig [3] 



be zero for any value of n xp , so that parametric oscillation 
can never take place in this case. The mean field shifts 
push in fact the signal/idler modes out of resonance be- 
fore the parametric coupling can overcome the damping 
rate 7 si . 



E. Laser intensity at threshold 

In the previous section we have determined the value 
of the pump mode population n xp at the threshold for 
parametric oscillation. The value of the corresponding 
laser intensity is then obtained by using Care has 
to be paid to the fact that single-mode instabilities may 
make some branches of the bistability loop dynamically 
unstable and therefore not reachable in an actual exper- 
iment. Again, this feature is typical^ of a x OPO and 
is absent in x*- 2 ' ones, where the relation between the in- 
cident intensity and the pump mode population in the 
pump-only state is a linear one and no instability other 
than the parametric one is possible^. 

Let us start from the 7 P = 7 S = 7,; = 7 case. The 
predictions for the value of the laser intensity at the 
parametric threshold are summarized in FigO where the 
contour plot of the threshold laser intensity is shown as a 
function of the detuning e s , — e p between the signal/idler 
mode frequencies and the pump mode frequency, and the 
detuning uj p — e p of the pump laser from the pump mode 
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FIG. 6: Pump-only characteristic curves for different val- 
ues of the detunings. Instabilities with respect to a para- 
metric oscillation process at a given e s i are indicated as 
red dashed, pump-only instabilities are indicated as black 
dotted. Arbitrary units for I p are the same in all panels. 
The labels of the plots correspond to the different regions in 
Fig[5] (A) parametric oscillation in the optical limiter regime 
— e P = —1-57, Lo p — e p = 0); ; (B) bistability regime: para- 
metric oscillation threshold above the pump-only bistability 
region (e s j — e p — —2*y, lu p — e p = O.97), the inset shows the 
bistability region more in detail; (C) bistability regime: para- 
metric oscillation threshold at pump-only bistability region 



= -O.67, 



e p = I.57); (D) bistability regime: 



parametric oscillation threshold not reachable with upward 
ramp in laser intensity (e s j — e p = — O.27, lu p — e p — I.57); 
(E) bistability regime: parametric oscillation threshold before 
pump-only instability point (e S j — e p = O.47, u) p — e p = I.57). 



frequency. Throughout all the present discussion, the 
laser intensity is assumed to be slowly but monotoni- 
cally increased from until the parametric threshold is 
reached. 

The lower right corner of this figure corresponds to the 
hatched region in FigU] where parametric oscillation can 
not take place because lo p is not sufficiently larger than 

The heavy horizontal line at ui p — e p — V07/2 separates 
the regions where the pump-only solution ^ respectively 
shows optical limiting (below the line) and optical bista- 
bility (above the line) . In the optical limiter case shown 
in FigJBJ^., the pump mode population n xp is a always 
a single valued function of the pump laser intensity I p . 
For a certain window in pump intensity, the (initially 
red-detuned) signal/idler frequency e S i is brought into 
resonance with the pump energy lu p by the mean-field 
shift, and the pump-only state becomes unstable with 
respect to parametric oscillation (red dashed line). Note 
that differently from the case of x*- 2 ^ OPOsi&, paramet- 
ric oscillation with x'' 3 '' media has an upper threshold as 
well: for too large pump laser intensities, the blue-shift 
of the signal/idler frequencies brings them off resonance 
and parametric oscillation can no longer take place. 

In the optical bistability case, the interplay between 
the pump-only hysteresis with the parametric oscillation 
leads to a variety of different behaviours (regions B-E). In 
order to fully understand these issues, it is useful to iden- 
tify the relative position of the pump-only and the para- 
metric instability regions on the n xp vs. I p curves which 
are plotted in FigJS] The different regions indicated in 
FigjS] correspond in fact to different arrangements of the 
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two instability regions. 

The simplest scenario is shown in FigJ5j3, where the 
signal/idler frequency e si is very red-detuned from the 
pump mode frequency e p . The pump mode population 
needed to bring the signal/idler modes on resonance is 
then much higher than the one needed to go through 
the pump-only hysteresis loop. In this case, paramet- 
ric oscillation occurs well above the bistability region so 
that pump-only bistability and parametric oscillation are 
effectively decoupled. The behaviour of parametric oscil- 
lation is completely analogous to the optical limiter case. 

For the parameters of Figj6]C, the pump only insta- 
bility still sets in before the parametric instability is 
reached, but the state of the upper branch where the 
system is expected to go, is parametrically unstable and 
OPO can start. This means that the laser intensity 
threshold for parametric oscillation coincides with the 
turning point of the hysteresis loop and in particular no 
longer depends on the signal/idler frequency e S i. For this 
reason, the contour lines shown in region (C) of Fig[5]are 
straight horizontal lines. 

FigEj3 shows a situation where the parametric oscilla- 
tion can not be reached by an upward ramp of the laser 
intensity. For increasing pump laser intensity, the sys- 
tem jumps to the upper branch of the hysteresis loop 
which is now parametrically stable in the region of inter- 
est, so that parametric oscillation does not start. Phys- 
ically, the jumps shown by the pump mode population 
at the switch-on point of the hysteresis loop is in fact 
large enough to make the signal/idler detuning to jump 
directly from one side to the other of the resonance. De- 
pending on the exact position of the parametric unstable 
region along the hysteresis curve, parametric oscillation 
can possibly be obtained by ramping the laser intensity 
down along the upper branch. Finally, Figj6fl describes 
the case when parametric instability sets in before the 
bistability saddle node bifurcation is reached. 

In Sec lIVI we shall see that the parametric instabilities 
shown in FiglBJ^-C lead to a stable OPO state. On the 
other hand, the situation is more complex for the case of 
Figj6j3, where it may happen that no stable parametri- 
cally oscillating state is available and the system even- 
tually ends up in the upper branch of the pump-only 
hysteresis loop. 



F. Quest for the lowest threshold 

In order to minimize the parametric threshold inten- 
sity, a careful choice of the detunings has to be performed: 
in this section, we will show that the mean-field shifts 
of the frequency modes make this optimization problem 
somehow more complex than a trivial question of "magic 
angle" . 

The optimal value of the detuning between pump fre- 
quency uj p and signal/idler frequency e s j is given by (fT3|) . 
In order to minimize the value of the incident pump in- 
tensity at threshold, one has to simultaneously impose 



a resonance condition between the pump laser frequency 
and the renormalized pump mode frequency: 



■g\x P 



(15) 



The optimal pump and signal/idler mode detunings are 
then immediately obtained by combining this result with 
111 and Q: 
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(16) 
(17) 



The corresponding value of the threshold intensity is ob- 
tained by simply substituting into Eq. ^ and then using 
Eq. (2}. For equal radiative and polaritonic decay rates 
Irad = 7, one obtains: 



NtrN QW h 2 j 2 uj p 



8CpXpX s Xi 



fJ 



(18) 



In Fig[51 the location of the minimum is indicated by 
a star in the present J s ,p,i = 7 case, this point lies 
in the region (A) where the behaviour of the system is 
the simplest. The pump-only solution being of the opti- 
cal limiter type, no hysteresis effects take place nor any 
interplay between parametric emission and bistability . 
Remarkably, both and A° pt and have a weak de- 

pendence on the Hopfield coefficients: for X SiP> i = 1, 



they are equal to A 



Opt 

si,p 



-7/2, A°p* = 7/2, while they 



are approximately A°f P = —0.537, A° pt = 0.527, for 
the typical values for a semiconductor microcavity used 
in Fig|5] This result is a refinement of the concept of 



"magic angle" at which perfect resonance A s 



is 



satisfied: as already noted in Ref|9|, small, but finite de- 
tunings A s j iP and A p are useful in a CW experiment to 
compensate the blue-shift of the signal/idler modes for 
increasing pump mode population. 

The value (j 1 T[) of the optimal detuning can be trans- 
lated in the wavevector space using the results of 
scc lIII Cl the optimal detuning is indicated in FigGJi by 
the horizontal line. It is easy to see that this value can 
actually be achieved as soon as k p > 1.34/mi - , which 
corresponds for the cavity parameters of Fig. [T]to a pump 
angle larger than 10 degrees. This minimum pump an- 
gle depends on the damping rate 7: for a smaller 7, the 
curves of Figj3] are stretched in the y-direction so that 
the optimal detuning can be already obtained at smaller 
pump angles. 

The red-shaded area in Fig[5] indicates the values of 
detuning A s .; jP that are available for a pump wavevector 
equal to k p — 1.4 /mi -1 . As the pump frequency uj p can 
be chosen at will, no bound exist in the vertical direction 
and this area is bound only in the horizontal direction. 
For given values of lo p and k p , the parametric oscillation 
dynamics is initiated when the incident pump intensity 
starts exceeding the minimum value of the threshold on 
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the horizontal segment contained in the red-shaded area 
in FigEl 

A crucial role in the OPO operation is played by the 
damping rate. From Eq. (|18[) . one sees that the laser in- 
tensity at the optimal point is proportional to the square 
of the damping. Furthermore, the value of the damping 
affects the extent of the red shaded area of available fre- 
quencies: for fixed pump angle, the border of this area 
moves to the right upon increasing the damping rate and 
eventually no longer overlaps with the (A) and (C) re- 
gions which are most favourable for OPO operation (see 
Sec|lV|). 

G. New features of the general 7 S 7^ 74 7^ 7 P case 




[ Gsi £p ] / Yp 



FIG. 7: The same as Fig[5]for unequal damping rates 7 P = 
7s = 7i/5. The lowest contour line in the plot is at 1.05/,™™! 
where the minimum of the threshold intensity I™™ is attained 
at the point indicated by a star *. The difference between the 
contours is 0.11™™. 



In current experimental configurations, e.g. for a pump 
in the vicinity of the "magic angle" , the idler linewidth 
is often much larger than the signal and pump ones, i.e. 
7; S> 7 P — 7s- Although the general formalism intro- 
duced in the previous sections keeps holding its valid- 
ity, some of the physical conclusions of Sec lIIIFl have to 
be modified. Because of the increased damping rate of 
the idler modes, the pump mode population n xp required 
by (flD)) for the onset of the parametric oscillation corre- 
sponds to a blue-shift of the pump mode g\X p \ 2 n xp much 
larger than j p . With the frequency choice suggested by 
a naive application of the condition (| 1 5[) . the intensity 



value (fl8|) would correspond to the end-point of the up- 
per branch of the bistability curve. Unfortunately, this 
point can not be reached by the simple upward ramp of 
the pump laser intensity considered in the present paper, 
so that a more complete analysis is required which fully 
takes into account hysteresis effects. 

The results are shown in Fig [7] for 7 P = 7 S = 7i/5- 
Because of the high value of the pump mode blue shift 
at the onset of parametric oscillation, the (D) and (E) 
regions are shifted to large values of A p (not shown), 
direct contact between the regions (A) and (C) is lost 
and the gap is filled by the (B) region. The optimal 
choice of the detunings lies on the border between the 
regions (B) and (C): parametric oscillation starts on the 
upper branch of the pump-only hysteresis curve exactly 
at the landing point of the jump from the lower branch. 
As one can see on Figj7l the optimal values of the A"^ 
and A° pt detunings (measured in units of 7 P ) are here 
larger than in the previous 7 P . s ,i = 7 case. On the other 
hand, the threshold intensity is increased above the naive 
prediction (fT5|) by a moderate factor of the order of two. 



H. Quantitative discussion 

In many technological applications of optical paramet- 
ric oscillators, a value as low as possible for the threshold 
intensity can be a key advantage. In this respect, semi- 
conductor microcavities in the strong coupling regime 
are very promising systems thanks to the extremely high 
value of the nonlinear coupling constant <?, much higher 
than the one of OPOs based on passive x^ media (see 
discussion below Eq. ((4])). 

Using typical values fi.7 = 0.1 — 0.5 meV for the damp- 
ing rates and Nqw — 3 for the number of quantum 
wells inside the cavity, Eq. (fT5)) yields a value in the 
0.13 — 3.2kW/cm 2 range for the incident laser intensity at 
the parametric threshold, a value which is in rough agree- 
ment with experimental data of Refs. Given the 
scaling relation (fT!?|) , the threshold intensity of passive 
^•( 3 ) OPOs with a comparable quality factor is orders of 
magnitude higher. Unless planar cavities of much higher 
quality factor are developed to compensate for the much 
weaker nonlinearity of passive materials, semiconductor 
microcavities in the strong coupling regime appear to be 
most favourable systems in view of low-power OPO ap- 
plications. 

In order for this comparison to be fair and complete, 
it is important to extend it to the case of OPOs based on 
passive x^ 2 ' material a 16 ! 17 ' 18 . In this case, no mean-field 
shift of the mode frequencies occurs and the minimum 
value of the parametric threshold is attained under the 
resonance condition u> p — e s , and is equal to: 

\P\min= 7/(232), (19) 

where the second-order nonlinear coupling constant for 
a planar cavity of thickness d and filled of a medium of 
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linear dielectric constant en n is 

%2=CX<2) \/W- (20) 

At optimal pump detuning, the driving amplitude equals 
\Fp\ = 7 1 P 1/2; irrespective of the type of nonlinearity. 
Combining this result with eqns. ©,6]), (H2), (HU) and 
(j20f . one finds the ratio between the threshold laser in- 
tensities of an d x' 3 ' OPOs: 

r(2) ~ y(3) • ^ 

Using a typical value x^ — 10~ 9 esu for a large Kerr 
nonlinearity, the value x^ 2 ' = 4 x 10~ 8 esu of the widely 
used KTP crystal 2 !, and j/u p = 2x 1(T 4 as in a typical 
A/2 semiconductor microcavity, the ratio (f2Tj) turns out 
to be around 0.016. 

This argument concludes the verification of the 
widespread expectation that for comparable values of the 
quality factor, the threshold intensity for parametric os- 
cillation is orders of magnitude lower in semiconductor 
microcavities in the strong coupling regime than in OPOs 
based on passive x^ 2:3 ^ materials. 

IV. BIFURCATION TYPE AND NONLINEAR 
SOLUTION ABOVE THRESHOLD 

In the previous section, our attention has been fo- 
cussed on the behaviour of the system below the para- 



metric threshold and we have characterized the value of 
the threshold intensity as a function of the detunings. To 
complete the study, it is then necessary to investigate the 
nature of the threshold point and characterize whether 
the onset of parametric oscillation takes place in a con- 
tinuous or discontinuous way. Both kind of behaviours 
have been indeed observed in the experiments 5 - ' 15 ' 23 . Cor- 
respondingly, the theoretical analysis in the present sec- 
tion will show that a variety of different hysteresis effects 
can take place depending on the kind of bifurcation that 
occurs at the threshold. Again, our discussion here will 
be based on the assumption that k S i are given quanti- 
ties. A complete discussion of the selection problem is 
postponed to a forthcoming publication^. 



To make the analysis the simplest, a three- mode ansatz 
of the form: 



V>Lp(k, t) = S k.k.e-*"'' + P S kMp e-* a > t + I S k , ki e~ iUit , 

(22) 

can be used&i 2 -, the signal/idler frequencies and wave 
vectors being related by Wi = 2ui p — lo s and = 2k p — k s . 
By projecting the wave equation (TT]) onto the three signal, 
pump, and idler modes, the following equations of motion 
are found 



4p = 

dt 


e p 


■ 7 

- l 2 


i—b = 
dt 




• 7 
- % 2 


d - 

i-rl = 
dt 


(, 


.7 
_ *2 



P + gX; ( |P| + 2|5| + 2|/| J P + 2P*SI 



S + gX 2 s 



2\P\ Z + \S\ Z + 2\I\*) S + P Z I 



2uj p + uj s 



i + gx 



(2|P| 2 + 2|£| 2 



ir)I + P 2 S 



1 Q* 



(23) 
(24) 
(25) 



where the following shorthand notations have been in- 
troduced e p , s ,i = e(kp )S) i) and X p s i — X(k p s i ). Scaled 
quantities S = X s S , P = X p P, I = X { I and F p = 
X p F p have been also defined. Imposing the stationarity 
of the solution and the condition that to s is purely real 
gives a set of 7 real equations (3 complex ones, plus 1 
real equation) which has to be solved for a total of 8 real 
quantities: the three amplitudes S, P, I and the (com- 
plex) parametric oscillation frequency u> s . The extra de- 
gree of freedom which is left undetermined corresponds 
to the U(l) signal/idler phase symmetry which is spon- 
taneously broken above the threshold 3 ^. 



FigOshows the behaviour of the pump n xp = |P| 2 and 
signal n xs — \S\ 2 mode occupations as a function of the 
incident pump intensity I p for different choices of pump 
laser to p — e p and signal/idler detuning e s , — e p . These 
plots exemplify the system behaviour in the most sig- 
nificant among the regimes studied in FigJHl Full lines 
indicate stable regions, dashed lines are the unstable 
ones^. Correspondingly, a numerical integration of the 
time-dependent equation of motion (|23H25p has been per- 
formed for a laser intensity which is continuously swept 
up and down through the parametric threshold. The re- 
sulting time-dependence of n xp and n xs is shown in Fig[9] 
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FIG. 9: Time evolution of the pump laser intensity (upper 
panels), pump intensity (central panels) and signal intensity 
(lower panels) for the same detuning parameters as in the 
panels A, C and E2 of Figg] 



I [arb. units] 



I [arb. units] 



FIG. 8: Pump and signal intensity as a function of the pump 
power in the different regimes. Black heavy full (red thin 
dashed) lines refer to the (un)stable solutions. 7s, p .; = 7 are 
taken here, but the results are qualitatively robust to a change 
of the 7's. Arbitrary units for I p are the same in all panels. 
(A) Optical limiter regime (ui p — e p = 0, e s ; — e p = — 1.57); 
(C) Optical bistability regime where parametric and pump 
only threshold coincide (ui p — e p = I.57, e s j — e p = — O.57); 
(El) Optical bistability regime where the parametric thresh- 
old precedes the pump only instability. The parametrically 
oscillating solution has a very small stable part (uj p — e p — 
I.57, e a i — e p — O.47); (E2) Same as El, with a different set of 
parameters such that parametric oscillation is here possible 
for a wide range of pump intensities. (u p — e p — 37, e S i — e p — 
27). The inset shows a magnification of the pump population 
in the OPO regime. 



for the most significant cases. 



A. Region (A): Optical limiter 

In the optical limiter case of FigJ5J\, both the pump 
and the signal mode occupations are continuous functions 
of the pump intensity. The transition is analogous to a 
second-order phase transition: the signal intensity is zero 
below and at the threshold and increases smoothly as a 
function of the pump power. In the language of nonlin- 
ear dynamics, this corresponds to a so-called supercritical 
Hopf bifurcation^. The corresponding time evolution is 
shown in the plots in the left column of Fig|5J both n xp 
and n xs have a smooth evolution in time which is imme- 
diately understood by following the curves of Fig(5]A.. The 
kinks correspond to the points where parametric emission 
switches on and off. 



It is interesting to compare this behaviour with the 
one of the x (2) OPO in the A Ai < 1 regime discussed 
in Ref. 18. As one can see in Figs J5] and the optical 
limiter regime corresponds to e si — e p < and uj p < e p . 
In both cases, the Hopf bifurcation is supercritical, and 
the populations are continuous functions of the pump in- 
tensity across the threshold. The behaviour well above 
the threshold is however completely different. In the x 
case there is an upper threshold as well, so that paramet- 
ric oscillation disappears for too large a pump intensity 
(not shown in the time-dependent plots). In the 
case, the parametric oscillation takes place instead for 
all values of the laser intensity above the threshold. As 
shown in Refill, for very high values of the incident in- 
tensity it becomes however unstable towards self-pulsing 
and chaotic behaviour. 

The behaviour of the system for the parameters of 
FigJBjB is completely analogous to the optical limiter case: 
the pump only bistability and the parametric oscillation 
indeed take place in an independent way. In the OPO re- 
gion, the behaviour of the pump and signal populations 
as a function of the incident intensity is therefore closely 
analogous to the one shown in FigJ5]A. 



B. Region (C): Optical bistability 

The physics turns out to be much richer whenever 
parametric oscillation and pump-only bistability take 
place in the same range of intensities. In the case shown 
in FigJSTj, the pump-only solution loses stability at the 
pump-only saddle node bifurcation. As the upper branch 
of the pump-only hysteresis loop is parametrically un- 
stable, the parametric oscillation sets in. As shown 
in Fig[5p, the solution connecting the lower and upper 
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threshold for OPO can be a complicate (multivalued) 
function^: typically, there are two stable branches (in- 
dicated with a and /3), which not always can be reached 
in a continuous way by means of a simple upward ramp 
of the pump intensity. To determine which branch is 
actually selected, the dynamics of the system has to 
be considered (central column of FigfS]). We have nu- 
merically found that the system jumps to the a branch 
as soon as the pump intensity exceeds the pump-only 
turning point. This a branch is then followed during 
both the upward and the following downward ramps un- 
til the saddle-node instability at the end of the branch is 
reached. At this point the system has to jump to another 
solution: our numerical simulations have shown that the 
/3-branch is dynamically selected, where parametric os- 
cillation still takes place with an even higher amplitude. 
Finally, when the saddle-node instability at the end of j3 
branch is reached, the system has no choice but to jump 
back to the lower branch of the pump-only solution where 
parametric emission is no longer present. 

It is important to stress that this analysis is based on 
a three-mode approximation: although this is certainly 
a valid description of a three-mode cavity, it may be not 
representative of all that can happen in a many-mode 
system such as a planar microcavity, where the a branch 
is often Eckhaus-unstable against changes in the signal 
wavevector. On the other hand, the branch turns out 
to be generally much more stable. A more complete dis- 
cussion of these issues will be presented in a forthcoming 
publication^. 



C. Regions (E1,E2): Optical bistability 

In the E region, the parametric instability occurs be- 
fore the pump-only one, and the corresponding Hopf bi- 
furcation is generally of the subcritical type^. Two sub- 
cases are to be distinguished. 

For e si ^ e P (FigEfll), although the instability is of 
the parametric kind, no stable parametrically oscillating 
solution exists for any pump intensity above the thresh- 
old. The parametric threshold is in fact very close to 
the pump-only threshold, and only a very small part of 
the OPO solution is stable (circle in Fig[8|). As this sta- 
ble region corresponds to intensity values in between the 
upper and lower turning points of the pump-only bista- 
bility loop, parametric oscillation can not be reached by 
any continuous intensity ramp, no matter its direction. 

For e S i ^> e p (Fig[8j32), the parametric threshold is in- 
stead sufficiently lower than the pump-only one for a sta- 
ble OPO state to exist and to be reachable by means of an 
upward intensity ramp: a stable parametrically oscillat- 
ing solution exists in fact for laser intensities extending 
from well below to well above the parametric instabil- 
ity threshold. However, as the bifurcation at the lower 
threshold is of the subcritical Hopf type, parametric oscil- 
lation sets in in a discontinuous way for an upward ramp 
in laser intensity. A time-dependent calculation is then 



needed to ensure that the system actually jumps from 
the lower branch of the pump-only hysteresis loop to the 
parametrically emitting solution rather than to the upper 
branch of the pump-only bistability loop. The results are 
shown in the right column of Fig[9] the switch-on of the 
OPO emission during the upward ramp is discontinuous, 
as well as the switch-off during the downward ramp. This 
latter corresponds to a saddle-node instability at a pump 
intensity slightly lower than the one of the subcritical 
Hopf instability. A new kind of hysteresis loop is there- 
fore present: parametric emission gives in fact a positive 
feedback to the pump mode population and two solutions 
(a pump-only one and a parametrically emitting one) are 
possible in a range of pump intensity values. The main 
difference with respect to the standard pump-only hys- 
teresis loop is that the higher turning point is here at a 
Hopf bifurcation rather than at a saddle-node one. 

Remarkably, this phenomenology can be related to an 
analogous one shown by a x^ OPO in the Ao Ai > 1 
regime of Reffl8l. Indeed, lo p > e p and e S i ^> e p in our 
(E2) region. The qualitative shape of the parametrically- 
induced hysteresis loop is indeed similar, with the main 
difference of the hysteresis loop having a here a finite 
size also in the n xp vs. I p plot and not only in the P 
vs. I p one. A qualitative analogy with the x OPO 
can be found in the (C) case as well: in addition to the 
topological similarity, the pump mode population is a 
very flat function of I p along the (3 branch, and the phase 
of the pump mode amplitude P in the (3 branch differs 
from the one in the lower branch of the bistability loop 
in a way very similar to the phase hysteresis shown in 
the A Ai < 1 case of RefUl. 



D. Considerations on quantum fluctuations 

All the discussion so far has considered the polaritonic 
field as a classical one, and therefore has neglected its 
fluctuations around the mean-field value. Before con- 
cluding, it is interesting to shortly address the behaviour 
of the quantum fluctuations in the different cases. The 
fluctuations around the pump-only solution below the 
threshold are mostly determined by the nature of the 
instability at the threshold point, i.e. whether this is a 
single-mode or a parametric one. The physics of the fluc- 
tuations around the three-mode solution |22|) above the 
threshold is instead more comple:»2&, and here we shall 
limit ourselves to a few, very general remarks. 

In regions (A) and (B), the onset of parametric os- 
cillation closely resembles a second-order phase transi- 
tion: the signal, idler and pump mode populations have 
a continuous dependence on the pump laser intensity. 
The overall behaviour as a function of the pump laser 
intensity is qualitatively identical to the one discussed in 
Refl35l as a function of the pump laser frequency: as the 
threshold point is approached from below, the magnitude 
of the quantum fluctuations of the signal and idler beam 
monotonically grows and eventually becomes very large 
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in the vicinity of the threshold where an eigenvalue of the 
stability matrix © goes to zero. The fluctuations being 
due to parametric creation of signal-idler polariton pairs, 
the signal and idler beams show significant quantum cor- 
relatio n 36 i 37 ' 38 . Above the threshold, the signal and idler 
fields have a finite mean-field amplitude which continu- 
ously starts from zero. Quantum fluctuations around this 
three-mode mean-field solution have a more complex be- 
haviour: a quite general fact is that the importance of 
the fluctuations is most important close to the threshold 
point, and then quickly decreases as one moves far from 
the threshold^. 

In the (E) cases, the behaviour is almost the same 
in the region below the threshold: the instability hav- 
ing a parametric nature, the quantum fluctuations (as 
well as the quantum correlations) in the signal and idler 
modes grow as the threshold is approached and become 
strongest in the close vicinity of the threshold point. On 
the other hand the behaviour above the threshold point 
is dramatically different: the onset of the parametric os- 
cillation (provided it really starts, as in case E2) is dis- 
continuous, and a completely different solution branch 
is selected (Fig|Hpl2). Furthermore, the landing point 
on the new branch is not necessarily in the vicinity of 
the end-point of the branch, so fluctuations are generally 
moderate. Yet, their magnitude becomes again large as 
one approaches the end-point of the branch where one 
eigenvalue of the stability matrix around the three-mode 
solution ([22]) tends to zero. 

In the (C) case, the behaviour is very different already 
below threshold: as the instability at the end-point of the 
branch has a single-mode nature, the quantum fluctua- 
tions in the signal and idler modes remain moderate also 
in the vicinity of the threshold point, while the pump 
mode ones grow very large as typical of optical bistable 
systems 3 ^. 

V. CONCLUSIONS AND OUTLOOK 

In this paper we have given a systematic classification 
of the behaviour of a triply resonant optical parametric 
oscillator based on a semiconductor microcavity in the 
strong coupling regime. Because of the nature of 
the collisional excitonic nonlincarity, the interplay of op- 



tical bistability and optical parametric oscillation makes 
the behaviour of these systems much richer than the one 
of standard OPOs based on passive x^ 2 ' nonlinear ma- 
terials, and a variety of different threshold behaviours 
can be found already within a simple three-mode the- 
ory. In agreement with recent experiments, depending 
on the specific value of the detunings, either a contin- 
uous switch-on or a discontinuous jump can be found 
for the behavior of the signal intensity at the paramet- 
ric threshold. The different behaviours have been classi- 
fied by means of the general theory of bifurcations, and 
a simple relation between the nature of the instability 
point and the behaviour of the quantum fluctuations at 
the threshold point has been pointed out. 

In order to minimize the threshold incident intensity, 
a rigorous and quantitative refinement of the "magic an- 
gle" criterion is provided which takes into account the 
mean-field shift of the modes due to interactions, as well 
as the possibility of hysteresis effects in the pump-only 
dynamics. A slight blue-detuning of the pump laser and 
a comparable red-detuning of the signal/idler modes with 
respect to the pump mode frequency turns out to be 
favourable in order to compensate for the mean-field shift 
of the mode frequencies. 

Generalization of the theory to the many-mode case is 
under way. In order to take fully into account the in- 
homogeneous spatial profile of the pump laser spot and 
the competition between parametric oscillation in differ- 
ent k s modes, techniques mutuated from the theory of 
pattern formation in nonlinear dynamical systems turn 
out to be of great utility. 
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